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Abstract 

We use time-dependent spin-density-functional theory to study dy- 
namical magnetic phenomena. First, we recall that the local-spin- 
density approximation (LSDA) fails to account correctly for mag- 
netic fluctuations in the paramagnetic state of iron and other itin- 
erant ferromagnets. Next, we construct a gradient-dependent den- 
sity functional that does not suffer from this problem of the LSDA. 
This functional is then used to derive, for the first time, the phe- 
nomenological Gilbert equation of micromagnetics directly from time- 
dependent density-functional theory. Limitations and extensions of 
Gilbert damping are discussed on this basis, and some comparisons 
with phenomenological theories and experiments are made. 
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The collective behaviour of spins of mobile, correlated electrons that gives 
rise to metallic magnetism is still only imperfectly understood. The problems 
arising in micromagnetics JI| or spintronics JJ, to mention but two subjects 
of current technological interest that involve itinerant magnetism, are ther- 
fore often dealt with within the framework of the frankly phenomenological 
Landau-Lifshitz-Gilbert equations, or semiphenomenological models whose 
foundation in material-specific microscopic theory remains incomplete. On 
the other hand, first-principles spin-density-functional theory (SDFT), which 
gives a more or less adequate account of the ground state, fails even to ad- 
dress directly the issue of time dependence of the magnetization. Clearly, 
under these circumstances it is worthwhile to inquire what physical insight 
may be gained by approaching the problem from the novel point of view 
afforded by time- dependent SDFT (TD-SDFT) [0, |J. Our aim here is to 
explore the potential benefits and difficulties of such enterprise. 

In the present paper we thus (i) explain why the popular local-spin-density 
approximation (LSDA), on its own, cannot account for the magnetic fluctu- 
ations at and above the Curie temperature, not even when combined with 
TD-SDFT; (ii) construct a new gradient-dependent density-functional that 
does not suffer from this problem of the LSDA; and (iii) use this functional 
to give a microscopic derivation of the phenomenological Gilbert equation 
of micromagnetics. This equation has recently been invoked for the phe- 
nomenological interpretation of experiments on magnetic nanolayers ||, but 
has never before been derived microscopically from density-functional theory. 

Let us first consider fluctuations in the paramagnetic phase. Among the 
low-lying excitations of a ferromagnet above the Curie temperature are spin 
fluctuations. Temporarily, in the neighbourhood of an atom, these fluctua- 
tions give rise to magnetic moments that average to zero on a sufficiently long 
time scale. A convenient way to describe this is by introducing the average 
magnetization 

m T (r) = — / dtm(r,t), (1) 
r Jo 

where m(r, t) is the space and time- dependent magnetization density ob- 
tained, in principle, from solving the time-dependent Schrodinger equation 
for the many-electron problem, and r is the averaging time. 

Consider now the behaviour of this average magnetization as a function 
of t. Right after a magnetic fluctuation has given rise to a magnetic moment 
on a given site, the average m T (r) will be different from zero. Moreover, this 
moment certainly persists until r is of the order of the hopping time to the 
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nearest neighbour, T^, which can be estimated by the inverse band width, 
h/W. On the other hand, for times much longer than the one characteristic 
for thermal fluctuations of the magnetic moment, Tfi uc , the average over all 
configurations of m(r, t) is zero. These two regimes, one with short-term lo- 
cal moments the other with no magnetic moments, exist in any ferromagnet 
above the Curie temperature. Depending on the values of band width and 
Curie temperature there may also be a third regime, characterized by aver- 
aging times satisfying h/W << r << Tfi uc . In this case electrons hop from 
site to site, but maintain sufficiently strong correlations to prevent immedi- 
ate destruction of the local magnetic moments by thermal fluctuations. This 
latter regime has been addressed by the disordered-local moment picture || 
in which the electronic hopping is treated with DFT, and the much slower 
behaviour of the local moments is dealt with by using statistical mechanics. 

We now ask whether the above, widely accepted physical picture can be 
recovered within TD-SDFT. On a formal level there is no reason to expect 
otherwise, since TD-SDFT is an exact transcription of the time-dependent 
many-body problem, subject only to quite weak v-representability conditions 
0, but the situation is different when one considers the approximations for 
the density-functional that are necessary to perform an actual calculation. 
The simplest time-dependent density functional is the adiabatic LSDA, or 
ALSDA This functional is a straightforward generalization of the con- 
ventional LSDA of ground-state DFT, from which it can be obtained by 
simply substituting the time-dependent densities n(r, t) and m(r, t) for the 
ground-state densities n(r) and m(r). 

The first objective of this paper is to point out that a treatment of mag- 
netic fluctuations within the ALSDA can never correctly account for the ob- 
served magnetic fluctuations. To understand why, we only have to consider 
the TD-SDFT equation of motion for the spin degrees of freedom. Neglecting 
spin currents, this equation takes the form [|J 



where B ext is the externally applied magnetic field, 7 the gyromagnetic ratio, 
and B^ c (r, t) is the exchange-correlation (xc) magnetic field of TD-SDFT. 
Within the ALSDA this field is calculated as 



<9m(r, t) 



7m(r, t) x [B«rt(r, t) + B xc (r, t)} , 



(2) 



dt 




T,t) 



5E LSDA 



[n(r), m(r)] 



(3) 



5m(r) 



n(r)— ♦n(r,t),m(r)— >m(r,t) 
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where E LSDA is the LSDA functional of static SDFT. Within the LSDA B xc 
is locally parallel to m and hence drops out of equation (0). In the absence 
of an externally applied magnetic field this equation then predicts a time- 
independent magnetization density m(r, t) = const. Within the LSDA the 
average magnetization m T (r) is thus independent of the averaging time, and 
the LSDA is seen to be unable to account for the three different regimes 
described above. The LSDA thus does not capture the dynamics of magnetic 
fluctuations. 

This does of course not imply that previous LSDA-based calculations of 
spin fluctuations |6|, |[, which have had considerable empirical success, are 
wrong. Such calcuations typically avoid the problem by bringing in concepts 
from outside of that framework, such as external constraining fields that are 
not recalculated selfconsistently, adiabaticity assumptions for the spin dy- 
namics, fitting to model Hamiltonians, identification of Kohn-Sham energies 
with excitation energies and the Kohn-Sham susceptibility with the many- 
body one, etc. However, it would clearly be desirable to have an approach 
that is exclusively based on TD-SDFT, without the need for additional ap- 
proximations and assumptions. This observation provides the motivation for 
the development, in the present paper, of a prototype density functional that 
is designed to avoid the problems of LSDA. 

To this end, note first that the problem is an intrinsic defficiency of the 
adiabatic LSDA, which persists even if the functional is formulated in a non- 
collinear way, i.e., in terms of the full magnetization vector m(r, t) instead 
of the more common variables n^(r,t) and ri|(r, t). In Ref. [|J we found that 
the simplest way to graft a nontrivial spin dynamics onto the functional is to 
include gradient terms. However, not any gradients will do the job. First of 
all, the gradients must involve all components of the magnetization vector. 
This discards all standard GGA-type functionals, which are explicitly formu- 
lated for the z-component only and can therefore not account for directional 
fluctuations of the magnetic moments. Second, the functional must satisfy 
the zero-torque theorem (ZTT), which states that the net torque exerted by 
B I(; on the system whole must vanish HI: 

y>r m(r,t) x B xc (r,t) = 0. (4) 

Guided by these considerations we now construct a simple model for a 
B xc functional that avoids the problems of the LSDA. Our starting point is 
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the following (still completely general) representation of B xc as 

B xc (r,t) = [ d 3 r' f dt'K[n,m](r,r',t,t')m(r',t'), (5) 

J J — OO 

where K[n, m](r, r', t, t') is an, as yet undetermined, tensorial kernel, charac- 
teristic for the spin-spin interactions. To make contact with the phenomeno- 
logical theories discussed below we now assume that K is a short-ranged 
isotropic function of the form K(\r — r'\, t — t'). This isotropy assumption is 
reasonable if B xc and m, and therefore also K, are interpreted as averages 
over intraatomic distances, as is the case in the phenomenological theories 
||. Due to the short-rangedness of K we can then expand the vector m(r', t') 
under the integral in Eq. (|5|) about the point r and the time t, 

m(r', t') = m(r, t) + 
(r - r)V ® m(r, t) + -[(r' - r) ■ V] 2 m(r, t) + (f - t)m(r, t) + . . . , (6) 

where we have kept first-order terms in the temporal variation, m(r, t) = 
<9m(r, t)/dt, and second-order terms in the spatial one. As will become 
apparent later, it is these orders that are required to make contact with 
phenomenological theories. Based on experience with similar expansions in 
static SDFT, nonlinear terms in m (which would correspond to cubic or 
higher-order terms in E xc ) are not included in this expansion. Ti xc then 
becomes 

B xc (r, t) = K (t)m(r, t) + K 2 (t)V 2 m(r, t) + A(t)m(r, t), (7) 

with 

rt roo A 

K {t)=4n dt' dxx 2 K{x,t-t r ) (8) 

J-oo JO 



-co 

ft 



rt poo A 

K 2 (t) = 27r dt' dxx 4 K(x,t -f) (9) 

J-oo Jo 
r-t poo A 

A(t)=4n dt' dxx 2 (t' -t)K(x,t -t'), (10) 

J-oo Jo 

and x = |r — r'|. Substituting this in the equation of motion (^) we obtain 

m(r,t) = -7m(r,t) x [B ext (r,t) + B xc (r,t)] = 
- 7 m(r, t) x [B e , t (r, t) + K (t)m(r, t) + i^(t)V 2 m(r, t) + A(*)m(r, t)] .(11) 
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Equation flTT|) constitutes a tensorial generalization of the Gilbert equation 
of micromagnetics fll], § . That equation is usually written as 

ct 

m(r, t) = -7m(r, t) x B e// (r, t) + ^m(r, t) x m(r, t), (12) 

where B e jj is a phenomenological effective field comprising external and 
exchange fields, M s is the saturation magnetization, and a is the so called 
Gilbert damping constant |5| . We see that our equation ( O ) reduces to an 
equation of the form (|T2|) if the kernels we assumed for generality to be tensors 
are taken to be scalars, and their intrinsic time-dependence is neglected. 

Interestingly, after the reduction of tensors to scalars the first two terms 
(containing Kq and K2) in B xc satisfy the ZTT identically On the other 
hand, the Gilbert damping term on its own is not guaranteed to satisfy 
the ZTT. Since the latter is an exact constraint that must be obeyed by any 
effective field to be used in the equation of motion (Q), this must be considered 
a defficiency of the simple form of the Gilbert damping or, equivalently, the 
linearization of the expression for B^r, t). Another limitation of Gilbert 
theory that becomes apparent from our derivation is that in general the 
coefficients in Eq. (|TT|) are time dependent, while those in Eq. (|12"D are not. 
This time-dependence can give rise to additional dynamics and damping 
that is not described by the Gilbert term. Below we will show that there is 
a simple model within which one can rationalize this lack of intrinsic time 
dependence of the coefficients. 

Carrying on, momentarily, with time-independent scalars Kq and K% in- 
stead of time-dependent tensors, we can write the B xc functional as 

B xc (r, t) = K m{r, t) + if 2 V 2 m(r, t) + Am(r, *), (13) 

where Kq, Ki, and A are numbers that characterize the spin-spin interaction 
kernel. The correspondence between the phenomenological and the micro- 
scopic equations, found above, shows one way in which these numbers can 
be obtained, since the parameters entering the Gilbert equation can be ex- 
tracted from experiment or simulations |1|: The first term in Eq. flT3| ) is of 
the form obtained in the LSDA, i.e., parallel to m. (This observation implies 
that one can, in principle, determine Kq from the LSDA.) Interestingly, the 
derivative terms in the expression (p~3|) imply that the resulting B xc is not 
parallel to m even when the coefficients are reduced from tensors to scalars. 
As a consequence they give, unlike the LSDA, rise to nontrivial spin dynamics 
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and spin fluctuations when substituted back into the equation of motion. To 
show why we expect the resulting spin dynamics to be essentially correct, in 
spite of the approximations made in the argument so far, let us, temporarily, 
neglect the damping term, and write B e //(r, t) = B e:ct (r, t) + i^ 2 V 2 m(r, t) 
(the term containing Kq does not contribute to the equations of motion in 
the scalar approximation). This is precisely of the form of the effective field 
entering the phenomenological Landau-Lifshitz equation , which is known 
to correctly account for the observed spin dynamics in ferromagnets. In this 
context K2 is usually replaced by the spin wave stiffness D, which is re- 
lated to our K2 by D = r yTiM s K2 and can be obtained from experiment or 
independent calculations (see, e.g., Ref. [11] and references therein). 

One way to completely fix the parameters in our functional is thus to 
determine Kq from the LSDA, K2 from the spin stiffness, and A from the 
Gilbert damping constant. Clearly, this empiricism diminishes the predictive 
power of the functional, and we will propose below a slightly less empirical 
strategy in its place. As it stands, the main virtue of the construction leading 
to Eqs. (|TT]) and (|I~3| ) is rather in its conceptual consequences: (i) it gives some 
degree of microscopic justification to the phenomenological approaches, and 
(ii) it shows what the term missing in the LSDA treatment of spin dynamics 
and spin fluctuations is, and how a simple prototype H xc functional that has 
this term might look like. 

The simple functional (|i~3|) is already sufficient to recover the two basic 
regimes for spin fluctuations above the Curie temperature, discussed in the 



introduction. To show this we substitute Eq. (|T3| ) into (|Tl~D and linearize by 
setting m(r, t) = M + n(r, t), where n describes small fluctuations about 
the equilibrium magnetization Mo- The resulting set of coupled differential 
equations for n is readily identified as describing damped oscillations with 
frequency 

a(q)= ^ w - k Jf + m, ,14) 

where huj(q) = h^(B + q 2 M s K 2 )/(l + a 2 ) = (g/2 B + Dq 2 )/(1 + a 2 ), B = 
\~B e xt\, and k(q) = 2au(q)/(l + a 2 ). Below the Curie temperature M is the 
dominating contribution to m(r, t) and to the average m T of Eq. ([!]). Above 
this temperature M = 0, and the r-averaged magnetic moment is entirely 
due to the damped spin fluctuations described by n(r, t). As a consequence 
of the damping, the average magnetic moment is zero for sufficiently large r. 
On a time scale shorter than that set by the damping constant A, however, 
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the contribution of n(r, t) to m r does not vanish, and one obtains a transient 
magnetic moment. This fundamental separation of time scales is not obtained 
from the LDA, within which A = and there is no damping. Furthermore, 
in the absence of the external field B eii the spin dynamics continues to be 
driven by the gradient term oc K 2 , whereas in the LDA K 2 = and there is 
no intrinsic spin dynamics at all (Q LDA = 0). 

We now proceed to give a somewhat more microscopic characterization 
of the coefficients Kq, K 2 and A. To this end we recall that according to 
Eqs. (g) to (0) these coefficients are all defined as certain integrals over the 
spin-spin interaction kernel K(\r — r'\,t — £'). A simple, but not unrealistic, 
model for the space and time dependence of this kernel is 

K(x,t)=Ke- x/x e- t/T } (15) 

which amounts to assuming that spin-spin interactions decay exponentially 
both in space and time. The amplitude K in Eq. fll^D is a measure for the 
strength of the spin-spin interaction, while A measures its spatial range and 
T its temporal range, i.e., its memory. With this model for the spin-spin 
kernel the integrals (§) to (|II]) can be evaluated analytically, and one finds 
Kq = 8ttKTX 3 , K 2 = 48vrA^TA 5 = 6X 2 K 0} and A = -8ttKX 3 T 2 = -TK . 
These relations have several interesting consequences. First of all, they re- 
place the coefficients Kq, K 2 and A, which had been introduced above in a 
purely mathematical way, by the physically meaningful quantities K, A and 
T. Second, although the kernel K(x, t) in Eq. QT5D has an explicit time de- 
pendence, the coefficients calculated with it depend only on the memory time 
T, and not on t itself. Within this model we can thus justify the replacement 
of time-dependent coefficients by static ones, made by the phenomenological 
theories [cf. our discussion of the limitations of conventional Gilbert damp- 
ing, below Eq. (O)]. On the other hand, recent experiments report enhanced 
Gilbert damping in Fe-Au nanolayers ||. In layered systems the isotropic 
exponential model ( jT5|) for K is, of course, unrealistic, and its main conclu- 
sion, the time-independence of Kq, K 2 and A, does not hold. The microscopic 
equation flTTp then retains a time-dependence in its coefficients, and attempts 
to recast it in the form of the traditional Gilbert equation (with a time inde- 
pendent A oc a) must result in the appearance of additional damping terms. 
Interestingly, additional damping has indeed been found necessary in various 
proposals || for explaining the experimental results. 

Returning now to the model fll5|) , we note that the expressions obtained 
from it for Kq, K 2 , amd A can be used, together with the above connection 
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of these coefficients with the phenomenological parameters in the Landau- 
Lifshitz-Gilbert equations, to deduce relations between the latter and the 
microscopic parameters K, A and T. An interesting example is 



TD / D 

~ V Gha ~ V 6jhM s K ' ^ ^ 

which shows how the microscopic quantities A and T (characterizing the spa- 
tial range and memory of the spin-spin interaction kernel) are related to the 
macroscopic parameters a (Gilbert damping), M s (saturation magnetization) 
and D (spin stiffness). To test the consistency of these relations we now con- 
nect them to another semiphenomenological approach, namely Stoner theory. 
Recall that in that theory there are no gradient-dependent terms, since it can 
be interpreted as a linearized LSDA. If we use it to evaluate A we thus ex- 
pect that the gradient- dependent terms in our functional disappear. This 
is indeed the case: In Stoner theory the proportionality factor between the 
effective magnetic field and the local magnetization is just 1/2/Iq times the 
Stoner parameter I. Hence K « 7/(2//q). Plugging this into the preceeding 
equation and using for I, D, and M s the experimental values for iron, we 
find that A = 0.18 x 10~ 10 m. This value, which is much shorter than the 
lattice constant in iron (2.7 x 10~ 10 m), shows that Stoner theory is indeed a 
local theory, in which spin-spin interactions decay very rapidly. It also shows, 
through the relation K 2 = 6X 2 K , that the gradient-dependent terms in our 
functional are strongly suppressed when one uses Stoner theory to evaluate 
them, as expected. 

Finally, we point out that the time average of our expression for B^, 
Eq. (|T5|), can also be interpreted as a magnetic equation of state fll2"| , and used 



to extract the equilibrium magnetization. Naturally, in the paramagnetic 
state the time- dependence of T5 xc (r,t) is due to fluctuations only, and the 
long-time average of B xc (r,t) is zero. In the ferromagnetic state, on the 
other hand, this average remains finite. Upon including a cubic term in m 
in the expansion and calculating the average according to Eq. (|l|), we indeed 
recover the usual equations of spin-fluctuation theory []12| |, determining the 



equilibrium magnetization and the Curie temperature. Here, however, they 
are obtained by invoking the ergodic theorem, characterizing equilibrium via 
temporal averaging instead of by thermal averages. 

The main results of this work are as follows: We have shown that the 
LSDA fails to account for spin fluctuations near the Curie temperature, even 
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when it is made time-dependent (ALSDA) and fully noncollinear. We next 
identified certain gradient-dependent terms as the key missing ingredient in 
the ALSDA, and constructed a very simple prototype functional that contains 
these gradients. Although simple, this functional gives rise to a still quite 
general equation of motion for the spin degrees of freedom, which can be 
identified as a microscopic version of the equations of motion of Landau- 
Lifshitz and of Gilbert. This identification provides microscopic support 
for these phenomenological theories and opens up a way for the empirical 
determination of the parameters in the functional. However, it also brings to 
light some shortcomings of the phenomenological approaches, which may be 
relevant for recent experiments. In summary, we have presented evidence that 
time-dependent SDFT may be a novel and useful alternative to constrained 
SDFT as a first-principles description of itinerant magnetism in metals. 
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